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Abstract 

The Wei-Norman technique allows to express the solution of a sys- 
tem of linear non-autonomous differential equations in terms of prod- 
uct of exponentials. In particular it enables to find a time-ordered 
product of exponentials by solving a set of nonlinear differential equa- 
tions. The method has numerous theoretical and computational ad- 
vantages, in particular in optimal control theory. We show that in 
the unitary case, i.e. when the solution of the linear system is given 
by a unitary evolution operator, the nonlinear system can be by an 
appropriate choice of ordering reduced to a hierarchy of matrix Ric- 
cati equations. Our findings have a particular significance in quantum 
control theory since pure quantum evolution is unitary. 



1 Introduction 



In two papers [T] and [2j written nearly fifty years ago, Wei and Norman 
developed a method for solving systems of linear differential equations with 
variable coefficients based on Lie-group techniques. As it is well known, 
solutions of non-autonomous linear systems of differential equations can be 
expressed in terms of time ordered exponentials. From this point of view 
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one can treat Wei-Norman formulae as a way to calculate such exponen- 
tials. More generally they provide local coordinates on a smooth manifold 
on which acts a (finite-dimensional) Lie group of transformations in terms of 
exponential map from the Lie algebra. 

Since then the Wei-Norman method has found numerous applications in 
control and system theory, see, e.g. [5], [1], [S], [S], [Z], [H], [H], [ID], as well 
as a basis for numerical approximation methods of integrating linear non- 
autonomous systems [11], [12], [13] . 

In recently rapidly developing theory of quantum optimal control, (see 
e.g. [IS], [16], [H], [18], [19]) of particular importance is the case when the 
underlying Lie algebra is su{N), since the pure quantum evolution preserves 
the norm of a state vector. We show that in this case the Wei-Norman 
method leads to a hierarchy of matrix Riccati equations, provided that a 
proper ordering of generators (basis Lie algebra elements) is chosen. Our 
findings hinge substantially on the structure of the su{N) algebra and its 
complexification sl{N,C), in particular on the properties of commutative 
ideals of the latter. 

In the next section we present briefly the basics of the Wei-Norman 
method for general complex semisimple Lie group. In Section [3] we sum- 
marize properties of the sl{N, C) relevant for our reasoning, in particular 
the structure of its abelian ideals. Results of this section are translated in 
Sections m and El to properties of the adjoint endomorphism and exponential 
function on sl{N,C). The final result for the Wei- Norman formulae in the 
case of 5l{N, C) is given in Section [61 The corresponding results for 5u{N) is 
easy to recover by restriction to anti-Hermitian generators. In Section [7] we 
give explicit formulae for low dimensions. 

2 General Wei-Norman method 

Let G be a n- dimensional Lie group and g - its Lie algebra. We assume in 
the following that g is complex and simple. Let also M 9 t M{t) G g be a 
curve in g and K{t) - a curve in G given by the differential equation 



In g we choose some basis X^, = 1, . . . , n in which M{t) takes the form 



d 
di 



K{t) = M{t)K{t), K(0) = /. 




n 




(2) 



k=l 
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and look for the solution K{t) in the form 

n 

K{t) = l[exp{uk{t)Xk), (3) 

k=l 

involving n unknown functions Uk- Differentiating we get 

n l — l n 

K' = "^u'lYl exp{ukXk)Xi Y[ exp{ukXk) 

1=1 k=l k=l 
n l—l 1 n 

= ]^exp(MfcXfc)X« Y\_ exp(-iifcXfc) ]^exp(MfcXfc) 

1=1 k=l k=l-l k=l 

n l~l n 

= ^u'lYl Adexp(nfeXfe) ■ XiK = '^u'iYl exp(iifcadxj ■ Xi K, (4) 

1=1 k=l 1=1 k<l 

where by ' we denoted the differentiation with respect to t and Ad is the 
adjoint action of G on g, 

Mg-X:=gXg-\ geG, X e Q. (5) 

In the last equality we used used Adexp(fex) = exp(6adx), where adx = [X, ■] 
is the adjoint action of g on itself. 
Comparing dl]) i 01]) we obtain 

n 

M{t) = ^ u[Y[exp{uka.dx,) ■ Xi. (6) 

1=1 k<l 

Both sides of (jH]) are elements of q and expanding both in the basis X^ we 
arrive at a system of coupled differential equations for Uk in terms of a/, 
k,l = 1, . . . ,n. 

It is worthwhile to rewrite (|6]) in more compact form. Let us denote 

A(') = JJexpKadxJ, (7) 

k<l 

and, consequently 

n 

A«-X, = ^4)X,. (8) 

i=i 

Now (El) can be written as 



n n n / n 



j=l 1=1 j=l j=l \l=l 



hence 

n 

1=1 

or in a compact form 

a = Au', (11) 

where A is an n x n matrix with elements Aji = A'-ji, i.e., its l-th column is 
equal to the l-th column of the matrix A^^\ c.f. ([S]), and a nad u are the 
vectors of the coefficients of M in ([2]) and the unknowns Uk- HA is invertible 
we obtain thus a system of (nonlinear) differential equations solved for the 
first derivatives 

u' = A-^a. (12) 



3 The structure of 5i{N^ C) algebra 

Let us start with a brief summary of the structure of simple Lie algebras [20] . 

Each complex simple Lie algebra g can be decomposed into the root 
spaces with respect to a chosen Cartan subalgebra i) (a maximal commutative 
subalgebra of q) 

S = f)©00a, (13) 

where the one- dimensional root spaces are defined as 

Qa:={Xeg: [H, X] = a{H)X MH G \)] (14) 

The linear forms a G f)* (the dual space to the algebra f)) are called roots 
and the element Xa spanning the subspace is called the root vector. By 
fixing a basis {Hi, . . . , if^} in P) we can identify the set of roots A with a set 
of r-dimensional vectors a = {a^, . . .a''), where = a{Hk). The number 
r = dim f) is called the rank of g. Among the roots we can find a set 11 of r 
positive simple roots, U = {ai, . . . a^}, such that for each a G A 

r 

a = '^miai, (15) 

i=l 

where either all mj are nonnegative (such roots a for a set of positive roots 
A+) or all nii are nonpositive (such roots a forming the set A_ of negative 
roots). There is one-to-one correspondence between positive and negative 
roots: for each positive root a the is a negative one —a. 
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The root spaces have the following property 

[0a,0/3] C Qa+p, (16) 

(for 13 = —a we have [ga,0_Q] C [)). 

In terms of positive and negative roots the decomposition f ll3l) can be 
rewritten as 

= n_©[i©n+, n±=00„. (17) 

aGA± 

The subalgebras 

b± = f)©n± (18) 

are called the Borel subalgebras relative to the Cartan subalgebra f). 

In the following we take G = SL{N,C) and g = 5i{N,C). In this case 
n = N"^ — 1. To construct a basis in q we define 

Ski = eke], (19) 

where e^, k = 1, . . . , N, is the standard basis in C^, ie., (e^); = Sm- The 
commutation relations for Ski read 

[Sij, Ski] = SkjSii — dilSkj. (20) 

For k 7^ I elements 5*^^ are the root vectors. The root corresponding to roots 
space spanned by Spq will be denoted by apg. From 5'^/ we construct a basis 
by the following renumbering 

= S„, ^^ iN + q-l)iN-q) fo,p<g^ (21) 

Xm = Sii-Si+,,i+i, m = ^N{N + l = l,...,N-l, (22) 
= V m = iV^- (^ + ^-f^-^^ p, forp<,. (23) 

One may choose roots apg corresponding to root vectors Spq for p < q to 
be the positive roots. Then the matrices (12T|) and (123|) generate maximal 
nilpotent subalgebras n+ and n_ respectively (conjugated with respect to 
the standard Hermitian structure on C^), whereas ( !22|) generate the Cartan 
subalgebra of sl{N, C). In this basis n+ is the set of strictly upper triangular 
matrices and n_ is the set of strictly lower triangular matrices. 

The set 11 of positive simple roots consist of elements of az,z+i for / = 
1,...,N-1. 



5 



In what follows we will show many useful features of this basis. First 
observe that according to f l^ . n+ := spanjXi, . . . , Xn(n-i)/2} and the order 
of root vectors spanning n+ corresponds to the following order of roots: 

>- Cik,i j < I or {j = I and i > k). (24) 

Analogously n_ := spa.ia{X]\f(^]\r^iy2, ■ ■ ■ , Xn2_i} and the order of root 
vectors spanning n„ corresponds to the following order of roots: 

cti,j >- <^k,i i > k or {i = k and j < I)- (25) 

Recall, that by the theorem of Lie [20], if L is a nilpotent Lie algebra of 
dimension r then there exist a sequence of ideals Ik fulfilling 

= Jo C /i C . . . C C = L (26) 

and a basis Yi, . . . , 1^ in L, such that 

4 :=span{Fi,...,n}. (27) 

The basis {Xi, . . . , Xiy(^iy_iy2} of n+ defined in (l2Ti) has this property. It 
follows from (ITB]) and the fact that aij + a^^i G A provided i = I oi j = 
k. In both cases aij + a^^i -< ctij, as can be easily checked. Thus Ik '■= 
span{Xi, . . . , Xk} is an ideal in n+ for any k = 1, . . . , ^N{N — 1). 

Direct consequence of Lie theorem is that in this basis the matrix of 
endomorphism adx^ : n+ — > n+ is strictly upper triangular |20]. More- 
over, according to fl25]) matrices generating n_ are transpositions of matrices 
generating n_ ordered in the reverse order. This implies that matrices of 
endomorphisms adxj. : n_ — )■ n_ for k = ^N{N + 1), . . . , A^^ — 1 are strictly 
lower triangular. 

We may decompose n+ in a particular direct sum of subspaces to modify 
slightly the original Wei-Norman expansion (|3]). To this end let us consider 
the family {Ofc}, k = 1, . . . , N — 1 of subspaces (the numbering of the root 
vectors Xi is the same as in (!?!]) ). 

k(k — 1) 

Ofc = span{Xi^^,Xi^+i . . . , Xi^^+N-k-i}, ik = N{k - 1) h 1. (28) 

In the standard matrix representation of sl{N, C) the subspace consists of 
matrices with only nonvanishing entries above the diagonal in the {N — k + 1)- 
th column. 

N-l 

We have dima^ = N — k and n+ = © afc. Each is an abelian subal- 

k=l 

gebra of b+ and Oi is an abelian ideal of b+. Let us also define 

N-l 

bfe = 0az. (29) 

l=k 
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Obviously C C bk-i- Moreover, each is a subalgebra of n+ and 
is an abelian ideal of b^, i.e., 

[bfc, bfc] C bfc, [ofc, Ofc] = 0, [ofc, bfc] C a^. (30) 

The subalgebras of n_ will be denoted by a^. Each is the Hermitian 
conjugate of and is generated by the basis elements ( 123|) in the following 
way: 

_ k(k + 1) 

ak = span{Xj^, Xj-^+i . . . , Xj^+N_k~i}, jk = N{N - k) ^ . (31) 

The subalgebras bk are defined in analogy to (l29l) and together wit they 
follow relations analogous to (l30i) . We have the following decomposition of 
g = sl(A^, C) into semidirect sum of 2N — 1 commuting subalgebras 

N-l N-l 

= [)©0a.©0H,. (32) 

i=i j=i 

4 Properties of the adjoint endomorphism 

Expression ([6j) consists of products of exponents of adx^ • In this section we 
present some fundamental properties of the adjoint endomorphisms corre- 
sponding to root vectors in sl(A^, C) which are crucial for usefulness of the 
Wei- Norman method in case of SL{N, C) group. 

Lemma 4.1 Let a be the root and G s((A^, C) he the corresponding root 
vector. Then: 

1. the image of (adx^)^ is equal to Qa O'^d 

ker ((adx„)^) = span {{Xi, Xjv2_i} \ . 

2. (adxj' = 0. 

Proof: First observe that statement 2 follows from 1. To prove the latter 
recall that [20] 

{Ho,, a + (3 = 0, 

iVa,/3^«+/3, a + /3eA, (33) 
0, otherwise. 
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where G f) corresponds to root a and A^^^^ is some constant. We have 

{-a{Hc,)X^, a + (3 = 0, 

Na,^Na,a+pX^+^+l3, a + a + /? G A, 
0, otherwise. 

(34) 

Since the condition a + a + /3GAis never fulfilled for sl{N,C), equation 
(IMj) implies that (adx^)^ sends any element of n+ U n_ into Qa- And since 
for H ei) 

(adxj' (H) = [X^, [X^,H]] = ~a{H)[X^,X^] = (35) 

the only element of basis fl21tl23p on which (adx„)^ takes the nonzero value 
is X_a, what ends the proof of statement 1. □ 
Observe that if we have two commuting matrices of given nilpotency order 
r then the sum of the matrices is also nilpotent of order r. It follows from 
the Jordan theorem [20j - the matrices can be expressed in Jordan form in 
the same basis and are block- diagonal with the same blocks of maximal size 
r — 1. Since and are commuting subalgebras, the lemma UT] yields: 

Corollary 4.2 If X e Uk or X then (adx)^ = 0. 

Lemma 4.3 Let a be the root and X^ G sl{N, C) be the corresponding root 
vector. In the basis defined by 1^711^) the matrices of adx„ for a G A+ are 
strictly upper triangular and for a G A_ are strictly lower triangular. 

Proof: Let Xa G n+. In the previous section we have mentioned the theorem 
of Lie, which has a direct consequence that the matrix of adx^:, : n+ — )• n+ is 
strictly upper triangular. So the sector of the matrix of endomorphism adx^ : 
g — )■ corresponding do n+ is strictly upper triangular. Since [n+, f)] C n+ 
the only nonzero matrix elements of adx^, in the sector corresponding to i) 
lie above diagonal. Finally we consider the action a.dx^{X 13) for Xjs G n_. 
The result is nonzero in two cases: Q: + /3 = Oora + /3GA. In the first case 
adx^iXp) G f) and the corresponding matrix element is above diagonal. In 
te second case it may be easily checked directly from the definition fl24H25l) . 
that a + (3 -< (3. This finishes the proof for X^ G n+. The proof for Xj^ G n_ 
is analogous. □ 

Lemma 4.4 Let Xq, G C st(X, C) or X^ G C 5i{N,C), where a is 
the corresponding root. The subalgebras ai, for I < k and the subalgebra 
bfc © f) © bfc o^e invariant subspaces of adx^ ■ 
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Proof: Let X^, G a^. We consider three cases: 

1. E ai, I < k. In this case G bi and Xjs G bi, because at C bfc C b/. 
Thus, Sidx^iX/s) = [Xa, X^] G bi. On the other hand G a/ and is 
an ideal in bi, so adx„(X/3) G a^. 

2. "Kg bfc © () © bfc. ^e have also Xa G b^ © () © b^ and the property 
adx„{Y) G bfc©f)©bfc follows from the fact that bfc©f)©bfc is subalgebra 
of g and this follows directly from the definition of bk and bk (see (l29l) ). 

3. Xi3 EcLi, I < k. According to definitions (1251) and fl^ the roots are 

a = {t,N -k + 1), i<N-k + l, (36) 

I3 = ^N-1 + I,j), N-l + l>j. (37) 

We have 

adx^iXp) = [Xa, Xp] ^ ^ a + /3 E A ^ i = j OT k = I. (38) 

Condition A; = / contradicts the assumption k > I. So the only possi- 
bility is i = j, what implies a + (3 = {N — l + l,N — k + l). Since k > I, 
we have N — 1 + 1>N — k + 1, so adx„ (-^z?) = Xa+p G H;. 

The same reasoning holds for X^ G Ofc. □ 
Since for a given k we have the following decomposition of g: 

g = ai © . . . © ak-i ® (bk®l)® b^) ©Ofc^i © ... © Hi, (39) 

one block 

the lemmas 14. ![ 14.31 and 14.41 yield 

Corollary 4.5 Let G C sl{N, C) (X^ G C sl(X, C) j be the root 
vector. In the basis [2T^2B\] the matrix of endomorphism adx„ is nilpotent of 
order 3, strictly upper triangular (lower triangular) and block diagonal with 
respect to decomposition (E^. 

5 Exponential function on commuting nilpo- 
tent subalgebras of 5i(N, C) 

The corollary 14.51 implies the following 
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Corollary 5.1 For X & or X the matrix o/exp(adx) is a quadratic 
polynomial in adx- Moreover in the basis l[EW^) the matrix 0/ exp(adx) 
is upper triangular (lower triangular for X G o^j and block diagonal with 
respect to decomposition (E^i- 

It is thus convenient to rewrite the equation ([6]) as a sum of terms corre- 
sponding to commuting subalgebras of decomposition (132|) . We define: 

Uk ■■= Yl exp(MjadxJ, Uk := JJ exp(MiadxJ, (40) 

where Jk = {ik, ■ ■ ■ + N — k — 1} is the index range defined in fl28l) . 
so for given k index i E Jk numerates generators of subalgebra a^, and 
Jk = {ifc, ■ ■ ■ Tik + N — k — 1} is the index range defined in (13T1) . so for given 
k index i E Jk numerates generators of subalgebra a^. For Z < — A; — 1 we 
also define: 

ik+l jk+l 

Vki ■■= JJexp(MiadxJ, Vm := JJ exp(MiadxJ, (41) 

i=ik i=jk 

where for given k index i runs over first / — 1 elements of Jk or Jk- We set 

Uo = Un = Vko = Vko = 1. 
Observe, that commutativity of subalgebras implies that: 

f^fc = JJ exp(MiadxJ = exp I ^ uMx, j = exp (^^(^(Y:^^j^u,x,y) ' (^2) 

so Uk equals exp(adx) for some X G ak (analogously Uk equals exp(adx) for 
some X G a^), so Uk and Uk have the triangularity and block diagonality 
properties implied by corollary 15.11 and the same statement holds for Vki and 
Vki- From lemma l4?T] we have also the following crucial property: 

Corollary 5.2 The terms quadratic in parameters Ui appear only when Uk 
or Vki o-cts on element of^k for the same index k and the terms quadratic in 
parameters Ui appear only when Uk or Vki icts on element of ak for the same 
index k. 

Let if, = ^N{N — 1) + 1 be the number of the first element of Cartan 
subalgebra according to ordering (12111231) . We define also: 

H := Y\_ exp(MjadxJ, Hi := Y\_exp{uiadx,), (43) 
«=«[) «=«[) 

where / G {1,...,A^ — 2} and Xi generate Cartan subalgebra f) (see fl22l) ). H 
and Hi are diagonal matrices and Hq = 1. 
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6 Wei— Norman method in a properly chosen 
basis 



The equation can be now rewritten in terms of operators defined in (HOl) . 
fl4ip and (1431) in the following simplified way: 

N-lk-l N-k-2 

M{t)=J2U^^- E Vk^<^^,X,,^.+ (44) 

k=l 1=0 i=0 
N-1 N-2 

j=l i=0 

N-1 N-lN-k+1 N-k-2 

+ UUrH-J2 n ■ E (46) 

j=l k=l l=N i=0 

On the other hand the matrix M{t) by definition ([2]) is a linear combina- 
tion of generators Xk with given coefficients ak{t). The first step in order to 
solve the system of differential equations is to solve linear algebraic equations 
for m'^. This is done by matrix inversion (Iimi2l) . We will show now, that 
the expression (1441146 p allows for separation of the full set of equations into 
sectors corresponding to decomposition fl32l) . 

The first term in (jUj) corresponding to A; = 1 is an element of Oi. It 
follows from corollary 15.11 and the fact that it is a sum of terms of the form 
adx^ for X G Oi and y G Oi. We move this term to the left hand side 
of equation (1441146 p . The rest of the sum, remaining on the right hand side 
has a common factor Ui, so we multiply both sides by Ui^. After this 
operations right hand side does not contain terms proportional to elements 
of ai, because is an composition of action of block diagonal operators with 
respect to decomposition (15^ for A; = 1 acting on generators of g starting 
from 02- The left hand side reads 

(N-1 \ / n N-1 \ 

- E ^^-1 = E «^(^)^^- - E ^^-1 ' 
i=i J \j=i i=i J 

(47) 

where we used the definition ([2]). Since there are no elements spanning ai on 
the right hand side, the first A^ — 1 components of fHTl) have to vanish. These 
N — 1 equations depend only on Ui for i < N (since Ui depends only on them) 
and for i < N (since V^i are upper triangular, what follows from corollary 
15.11) . Moreover Ui is a quadratic function of function Ui, so we end up with 
matrix Riccati equation for functions Ui for i < N. Once this system of 
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equations is solved and its solutions are substituted into equation (1441146 p . the 
terms corresponding to generators of Oi cancel, and we obtain reduced system 
of equations in smaller subspace of q. Next we apply the same procedure to 
the remaining sum on the right hand side of ( I44M61) . It starts now from 
term spanned by generators of 02. The same reasoning as for Oi applies 
and after moving the term in question to the left hand side and multiplying 
both sides by U2 we obtain a matrix Riccati equations for functions Ui which 
multiply generators of 02 in ([3]). We substitute the solutions of this system 
of equations and the main equation reduces again. We keep repeating the 
procedure described above until we obtain the solutions for all coefficient 
functions Ui corresponding to all a^. The solutions will always come from 
Riccati type equations. What remains from the sum (1441146 p on the right 
hand side of flTT]) after inserting the solutions and multiplying both sides by 
product of all U^^ is the following: 

N-2 N-lN-k+1 N-k-2 

(l.h.s.)=5^/7.<+,X,,+. + i^-E U ^ Y. (48) 

i=0 k=l l=N i=0 

First sum corresponds to Cartan subalgebra f). Since i) is commutative we 
have 

exp (adx) Y = Y, for X, F G 1), (49) 

so in fHHj) HiXi^^+i = Xj^+j, what follows from definition fH3|) . The second 
sum in fHHj) is a combination of generators of subalgebras only, because is 
an action of diagonal operator H and lower triangular operators Ui and Vki 
on those generators. Thus the functions -Uj corresponding to generators of () 
can be found by simple integration of solutions to previously found Riccati 
equations. After substitution of these solutions and multiplication of both 
sides of fHHl) by we are left with equation where on both sides there 
are only terms proportional to generators Xj spanning Hfc, because all other 
terms have canceled. In order to find the solutions for remaining functions Ui 
we proceed as follows. First we multiply both sides by which belongs on 
functions Ui corresponding do Oi only. Derivatives u[ on the right hand side 
appear only in the first term of the sum, because of the block diagonality of 
Vkj operators. Thus this terms separates form the rest and has to be equal to 
the action of U^^ on the left hand side. The fact that there are no generators 
of ttfc in the equation and coroUarv 15.21 imply that this set of equations will be 
linear in functions Ui. Once we solve it, the solutions may be substituted into 
equation, and terms proportional to generators of Oi will cancel. Then we 
multiply both sides by U2 and proceed in the same way obtaining the set of 
linear equations for functions Ui corresponding to 02. We keep repeating the 
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described procedure until we find linear equations for all remaining functions 

Ui. 

The above described procedure enables the conversion of highly nonlinear 
differential equation ([1]) into hierarchy of matrix Riccati equations and linear 
matrix differential equations. This procedure is effective in the sense, that 
it provides an algorithm that may be applied directly. The authors have 
written a program in Maple which performs this procedure for any given 
and tested its successful performance up to = 10. 

It is worth mentioning that the crucial ingredient for realizing the de- 
scribed algorithm in practice is the order of the generators (12111231) . Once 
this base is used for computations and the inverse in f|T2l) is successfully 
computed, the separation of the system of equations comes up automati- 
cally. For large N computation of the inverse is the part of the algorithm 
with the largest computational complexity, which scales with as A^!. It 
follows from the fact, that the matrix A in (|TT|) which is in principal of di- 
mension (A^^ ~ 1) X {N"^ — 1) has a special block upper triangular form with 
the largest block to invert of the size (A^ — 1) x (A^ — 1). The inversion 
has to be realized by Cramers rule which has the mentioned computational 
complexity. 

7 Examples 

7.1 5'L(2,C) 

For N = 2 the system f[T^ reduces to one Riccati equation, 

u[ = tti + 2a2^i ~ 0,3^1, 

and two equations for the remaining unknowns, 

u'2 = a2 — dsUi, 
u's = a3e2"^ 

which give U2 and u^ by simple integrations, once solutions of (150|) 

7.2 5'L(3,C) 

For A^ = 3 we obtain from fll2p : 

1. A system of two coupled Riccati equations 

u[ = tti + (2a5 — a^jui + aQU2 — o-s'^i " ajUiU2, (53) 
u'2 = a2 + a^ui + (04 + a^)u2 — a^uiU2 — ajul, (54) 



(50) 

(51) 
(52) 

are known. 
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which for further reference we will rewrite in the form 



u 



(1) — ^(1) + ^(i)"(i) + U(i)ufi)b(i), 



with 



and 



u 



(1) 



Ul 
U2 



«1 

a2 



(1) 



2q.5 — Q.4 Og 

as a4 + as 



(55) 
(56) 
(57) 



2. An equation for which reduces to a scalar Riccati equation upon 
substituting solutions of fl58l) 



""3 = («3 - 08^2) + (2a4 - as + as^i - a7U2) U3 + (a/tii - gq) ul, (58) 

3. Equations for the coefficients of the Cartan algebra generators which 
are solved by single integrations once solutions of (153|) and (158|) are 
known, 



^4 = a4 — a^us + a7 (miMs — U2) 
u'k = as — ttgWi — a7M2 



(59) 
(60) 



4. Equations for the coefficients of the generators of the second ("lower- 
triangular") nilpotent subalgebra 



u' 



-U4+2M5 



^8 ^ (ag + ar-usje" 
which are solved by simple consecutive integrations. 



(61) 
(62) 
(63) 



7.3 S'L(4,C) 

A similar structure emerges for = 4. The system of 15 equations f[T^ 
separates into: 

1. A system three coupled Riccati equations, 



u 



(1) 



^(1) + ^(1)^(1) + "(i)"(i)t>(i 



^(i)"(i) 



(64) 
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where 





Ml 








-ai5 


U(l) = 




' ^(1) - 


a2 




— ai4 




Us 




^3 




-ai3 



and 



(1) 



-ag + 2ag ai2 an 

fl4 — ttj -\- flg -|- CLg flxo 

0-5 fig (X7 -|- CLg 



2. A system of two coupled Riccati equations, 



u 



(2) 



where 





U4 




U(2) = 


U5 





^(2) + ^(2)^(2) 



04 - 015^2 
«5 - «15% 



U(2)U(2)b(2)5 



(65) 



(66) 



(67) 



'(2) 



-ai2 + auUi 
-«ii + aia^^i 



(6^ 



and 



(2) 



-ay + 2as — ag — ai4U2 + ai^ui 
ttQ - auUs 



aio — ^2ai3 — ai5 
ar + as - ag - aisUs + ai^Ui 

(69) 

Observe that once a solution of ( 164|) is known, the system (1671) is closed 
since ( 168|) and ( 169|) are given in terms of some known functions. 

3. A scalar Riccati equation, 

u'q = ae - ai2M5 + auUiu^ - auUs 

+ (2 ay — as — ^5011 + U4,ai2 — Usai3 + ai^uiu^ — auUiu^ + 'U2ai4)it6 
+ (-aio + W4aii - arsUiU4 + aisU2)ul, (70) 

with coefficients depending on solutions of and f E7|) . 

4. The remaining 9 equations for Ur,...,ui5 which are solved by sin- 
gle interactions of functions constructed from the initial coefficients 
ay, . . . , ai5, and solutions of (^^, (jSZl) and (170]) . 
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